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We show that the running of operators which mix under renormalization can be computed fully 
non-perturbatively as a product of continuum step scaling matrices. These step scaling matrices 
are obtained by taking the "ratio" of Z matrices computed at different energies in an RI-MOM 
type scheme for which twisted boundary conditions are an essential ingredient. Our method allows 
us to relax the bounds of the Rome- Southampton window. We also explain why such a method 
is important in view of the light quark physics program of the RBC-UKQCD collaborations. To 
illustrate our method, using n/ = 2 + 1 domain- wall fermions, we compute the non-perturbative 
running matrix of four-quark operators needed in if — > tttt decay and neutral kaon mixing. Our 
results are then compared to perturbation theory. 



Introduction 

Lattice QCD has now reached a stage of high precision 
in flavour physics and plays a crucial role in the quest 
for new physics. Precision phenomenology requires using 
a non-perturbative scheme to renormalize quantities 
obtained from lattice simulations, avoiding ill conver- 
gent and low order lattice perturbation theory. One 
possibility is to use the Schrodinger functional (SF) 
scheme, a theoretically appealing method which allows 
for a smooth connection between low energy - where the 
hadronic matrix elements are computed on the lattice 

- and very high energy, where the perturbative series 
is accurate [J H]. This connection is done through 
the use of step-scaling functions [3H5]. Unfortunately, 
the range of operators actually computed in the SF 
scheme is rather limited; perhaps this is due to the fact 
that a peculiar perturbation theory (different from the 
infinite volume one) is required. Another possibility 

- very popular in the lattice community - is to use a 
RI-MOM type scheme [6] : such a scheme is theoretically 
sound, relatively easy to implement and the connection 
to MS or any other perturbative scheme is done using 
continuum infinite volume perturbation theory. Thus 
one benefits from the recent multi-loop computations 
achieved by several groups (see for example [ZHS])- In 
an RI-MOM type scheme one numerically computes 
the off-shell amputated vertex function Go of the 
operator of interest O between external states with 
given momenta and fixed gauge gluonic configurations. 
We project this quantity onto its Dirac-colour structure 
and obtain the quantity Aq (the choice of projector 
is in general not unique and the precise definition of 
the scheme depends on this choice of projector). We 
take all the quark masses to be degenerate and equal 
to m. Then one requires that in the chiral limit the 
renormalized projected-amputated operator matches its 
tree-level value [B], i.e. if a is the lattice spacing, n is 
the number of quark fields, ip, in the operator and fi is 
the renormalization scale (which depends on the choice 



of external momenta), and F is the tree level value of 
Aq, one imposes 

Z^{fi,a) Z:,''^'^{p,a) liia Ao{fi,a,m) F . (1) 

In this way one obtains, non-perturbatively, the scheme 
and scale dependent renormalization factors Zq(^, a). 
The precise definition of the lattice scheme S depends 
on the details of the implementation. In the second step, 
one converts the result to a scheme more appropriate for 
phenomenological applications. In the case of an effec- 
tive low-energy theory one can express the Hamiltonian 
by a sum of local operators multiplied by some Wilson 
coefficients. One must find a scheme where both the 
renormalized matrix elements of these operators and the 
corresponding Wilson coefficients can be computed. Typ- 
ically, one matches the Z factors obtained in the lattice 
scheme at /i '-^ 2 — 3 GeV to MS. Equivalently, (if pertur- 
bation theory is accurate) one can divide the vertex func- 
tion computed in the lattice scheme by the corresponding 
running to obtain the renormalization group invariant Z 
factor and then convert to the desired scheme and scale. 
It is important to note that the running and the match- 
ing are usually computed only perturbatively (we discuss 
some recent implementations of a non-perturbative run- 
ning in the next section). 

In order to keep the discretisation effects under con- 
trol, and at the same time access a region where per- 
turbation theory can be applied, RI-MOM schemes re- 
quire the existence of the so-called Rome-Southampton 
window: ideally one would require Aqcd M 
In practice, this window can be quite narrow or even 
closed and most of the lattice computations are done with 
a/i ^ 1. In fact the size of the discretisation effects de- 
pends on the lattice action, and it is usually admitted 
that for an off-shell 0(a)-improved action it is sufficient 
to require (a^/ir)'^ ^ 1. Concerning the other side of 
the window the situation was greatly improved when it 
was realised that a different choice of kinematics, called 
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non-exceptional, strongly suppresses some infrared con- 
tributions to the vertex function [TOl[Tl]. Nevertheless, 
this window is still an important limitation, in particular 
when considering light pion masses where the physical 
volume - and thus the lattice spacing - has to be large. 
This could become a real issue when one tries to tackle 
the very challenging computation of if — > tttt decays be- 
cause one needs even larger physical volumes. A new 
step scaling method has recently been introduced [12] 
that allows to address the window problem for RI-MOM. 
This method has already been applied for multiplica- 
tively renormalized operators [131 114) . 

Strategy 

In this work we apply the method presented in [12] , gen- 
eralised to the case where several operators mix under 
renormalization. The main idea is to separate the scale 
and the lattice spacing(s) at which the bare matrix el- 
ements are computed from the ones at which the run- 
ning and the conversion to the perturbative scheme are 
performed. Let us consider the case of four-quark oper- 
ators mixing: following again [6J, we define a matrix of 
amputated, projected vertex functions = Pj{Goi}, 
where the Oi form a basis under renormalization and the 
Pi projects and traces onto the Dirac-colour structure of 
the corresponding operator Oi. If we denote by F the 
tree-level value of A, equation ([T]) becomes 



X Imi 



ik 



(2) 



where, for convenience, we use the vertex function of the 
axial current Aa to fix Z^. Equation ^ defines Zfj 
where the scheme S depends on the choice of kinematics 
and projectors of the four-quark vertex function and on 
the choice of quark wave function renormalization |13j . 
At finite lattice spacing a and for a given renormalization 
scale fj, we consider the matrix 

a)= lim \A\{fi,a,m) A~^{p.,a,m)] , (3) 

m— >0 

and we define the step scaling matrix by 

tT'^(/i, s/i) = lim Y,^ {fi, sii,a) = lim \Rs{ii,a)RZ^{sfi,a)] 

a— >0 a— 

(4) 

We also note that ^A(a) cancels out in the ratio since, 
to a very good approximation, it does not depend on the 
scale ^. One important point is that although the quanti- 
ties A, Z and R depend on the details of the computation 



this is not the case for the step scaling matrix which has 
well-defined continuum limit and depends only on the 
choice of renormalization scheme S (and on the number 
of flavours) . In this work we use a scheme which is called 
(7^, 7^)-scheme in ,13j : it involves non-exceptional kine- 
matics with a symmetric point. Such a choice greatly 
suppresses unwanted infrared effects. 
The strategy that we are proposing can be summarised 
in the following way: 

• Consider a set of simulations with a rather large 
volume of spatial extent Lq where physical pion 
masses can be simulated and with moderately 
large lattice spacings oq, but small enough for the 
Symanzik expansion to converge. Compute the 
bare matrix elements of interest {O^'^'^'^ (ao)) and 
renormalize them in a lattice scheme S at the 
low energy scale /io, i.e. compute (C^(/io)) = 
lim rzg(^o,ao)(C''"°(ao))]- Of course the scale 

ao-fO 

^0 should be such that the associated discretisa- 
tion errors are small but, compared to the Rome- 
Southampton window, we do not require the non- 
perturbative effects to be small. Instead one just 
has to ensure that the finite volume effects are neg- 
ligible, so the renormalization window becomes 



io < /"o < {t^/o-oY 



Iterate the following step, with i = 1, 2, . . . , n: con- 
sider a set of simulations, with a physical volume of 
space extent Li < Li^i and a set of lattice spacings 
Oj. With the requirement that, on each lattice. 



is satisfied, compute the step scaling matrix by 
evaluating eq. Q. 

• At a scale /Zn high enough to apply perturbation 
theory, multiply by the perturbative matching and 
running factors corresponding to the desired scale 
and scheme (typically MS at fj, = 2 or fj, — 3 GeV). 
In this volume the usual Rome-Southampton con- 
dition holds 



^^QCD 



< fil < (vr/fln)^ 



In summary, the general equation can be written as 



(0^^'=(Ai))=C^'=^'^(/^)x[/'^(/i,A*„)xa'^(/i„,/i„_i)xa'^(Mn-i,Ai„-2)x...xa*(Mi,Mo)x (O^'M) , (5) 

^ V ' ^ V ■' 

Fine lattices Coarse lattices 



where C^^^^l/i) represents the matrix of matching factors which converts the Z matrix computed in the 
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scheme S to the scheme MS, and /7'^(/z,//n) is the 
usual running matrix in the scheme S computed in 
perturbation theory. In the previous equation we have 
made explicit the fact that our method consists in 
re-expressing the running matrix, usually computed 
in perturbation theory, by a product of continuum 
non-perturbative step scaling matrices. 

One notices that in the first n— 1 steps, the lower bound 
of the Rome-Southampton window (Aqcd) is replaced by 
a more advantageous one {L~^). In other words: we do 
not need to be in the perturbative regime. Moreover a 
better control of the upper bound is achieved by tak- 
ing the continuum limit of the step scaling function (in 
particular there is no discretisation error of order ap/i^ 
in contrast to the "naive" RI-MOM implementation). A 
couple of remarks are in order: 

• Using twisted boundary conditions circumvents the 
Fourier mode constraints and allows us to fix the 
orientation of the momentum while changing its 
magnitude. As emphasised in jl2j . thanks to this 
property we can compute the vertex functions for 
an arbitrary number of points lying on the same 
scaling trajectory. The continuum limit of the ver- 
tex function as a function of is then properly de- 
fined and in particular we do not need to perturba- 
tively subtract any lattice artefact. The continuum 
extrapolation is also easier since, if we know the lat- 
tice spacings with sufficient precision, we can sim- 
ulate any arbitrary momentum. This is also useful 
because our method requires that each momentum 
p = jjLQ, Hi, . . . has to be common to different sets 
of lattices. 

• Since we impose periodic boundary conditions, one 
could worry about the use of perturbation theory 
in a small volume, where the space extent is of 
the order of (or smaller than) Aq^j-, [TH]. How- 
ever we claim that in a non-exceptional graph with 
hard external momenta, decoupling will ensure, for 
fjL ^ , that our computation is free from fi- 
nite volume effects and finite volume perturbation 
theory is not needed. While this might need fur- 
ther investigation, it is not relevant here since we 
consider only "infinite" volumes (of spatial extent 
much larger than AqJ,q). 

Before closing this section we wish to mention other 
works on non-perturbative running in RI-MOM. Taking 
the continuum limit of the ratio of Z factors at differ- 
ent energies in an RI-MOM scheme was first proposed 
in pi)| but the authors did not address how to match 
the momenta computed with different lattice spacings 
such that the lattice artefacts have an a? expansion. Zh- 
estkov [17] used fine tuning of /3 in the quenched ap- 
proximation to exactly match the Fourier modes but did 
not define a quantity which has a well-defined contin- 
uum limit. More recent work has looked at the ratios of 
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FIG. 1: Example of box diagram contributing to if" — if 
mixing in the Standard model. 



Zs [T4l[T8l[T9]. However, how to continuum extrapolate 
the distinct Fourier modes with distinct lattice artefacts 
must be addressed. Some cases have model input like 
perturbative subtraction of the lattice discretisation ef- 
fects, or rules of thumb such as use of sin(p) instead of 
p. Instead here we follow [T2] and implement twisted 
boundary conditions to keep the orientation of the mo- 
menta fixed and give the observables a smooth a? depen- 
dence, and allow for a full non-perturbative continuum 
extrapolation. 

Renormalization of kaon weak matrix elements. 

In this section we give the definitions of the kaon 
four-quark operators that we consider in this work. We 
refer the reader who would like to find more details 
about this part to the recent reviews [501 HI] ■ 

In the standard model, neutral kaon mixing is dom- 
inated by box diagrams like the one shown in fig- 
urejlj The non-perturbative contributions are given by 

{K lOyy^^^lif"), where Oyv+AA parity conserv- 

ing part of (s7^(l — 75) d)( 37^1(1 — 75) d). It is well known 
that this operator belongs to the (27,1) representation of 
5t/(3)L X S'J7(3)r and renormalizes multiplicatively. To 
study neutral kaon mixing beyond the standard model it 
is useful to introduce the so-called SUSY basis of As = 2 
operators. In this basis Of^*=^ is the standard model 
operator, Of^*^^, i > 1 are beyond the standard model 
(BSM) operators. Denoting by a and /3 the colour in- 
dices, one has 



(27,1) Of ^ 

0^=^ = {Mi-iM{M^-i^)dp). 

(Sa(l - 75)^^/3) (S/3(l - lb)da) , 



(6,6) 



qAs=2 



= (Sa(l-75)cfa)(s/3(l+75)c^/j): 
= (Sa(l - 75)0^/3) (5/3(1 +75)rfa) ■ 



The operators have been studied with various lattice for- 
mulations, see for example j 22H24] . As we wrote explic- 



itly in the previous equations, Of '^^^ and O^^^^ trans- 
form hke (6,6) under S'J7(3)l x SU{2>)yi and then mix 
together under renormalization 



O 



As=2 



belong to 



Likewise Of''~^ and 
and also mix together. Thus in 



a scheme which preserves chiral symmetry the five-by- 
five renormalization matrix is block diagonal: the only 
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non-zero Z factors can be divided in three subgroups: 
a single factor for the (27, 1) operator and two, two-by- 
two matrices for the BSM operators. In practice it is 
convenient to work in another basis where all the opera- 
tors are colour unmixed and wc consider only the parity 
even component of the four-quark operators. Using the 
notation F (g) T — > {sa^da) (s^Fc?^) we define the renor- 
malization basis by: 



(27,1) 






7m + 


7m75 € 


5 7m 75 


(8,8) < 




= 7m® 
= 1(8)1 


7m - 

- 75 


7m75 
<8 75, 


5 7m 75 


(6,6) . 




= I8)I 


+ 75 

S) cTm"^ 


8)75, 





The explicit relations between the two bases are given 
in the appendix. We denote by Z'^"^^ the (block diag- 
onal) renormalization matrix defined in the renormaliza- 
tion basis. 

It is interesting to note that the renormalization factors 
of some As = 1 operators which appear in iiT tttt de- 
cays can be obtained from those of the As = 2 operators 
mentioned above. At low energy in the A/ = 3/2 channel 
they are three operators that contribute: a (27, 1) called 
renormalizes multiplicatively and 
two (8, 8) which mix together: the electroweak penguins 

We give the explicit form of these opera- 
tors in the appendix. Denoting by Z^'^^^^, i, j = 1, 7, 8 
the corresponding renormalization factors, the relation 

to the Z,^"=2 reads 
''J 



^87 



'7As = 

^11 
^22 



= -2Zk 



^11 , 
5'As=l 
^78 
-3'As=l 



1 7As=2 

2 ^23 : 

^33 



In this work we give some results for the non-perturbative 
running of Q^^j^^ 5 and Q'f 7'^^''^^"^^^. A full computa- 
tion of _ftr — )■ TTTT decays requires also the renormalization 
of A/ — 1/2 four-quark operators which transform 
like (8, 1) under S'J7(3)l x SU^i)^. Because for these 
operators one has to compute disconnected diagrams we 
do not consider them in this work and leave them for the 
future. Although these operators are important in a full 
computation of if — >■ tttt amplitudes the main purpose 
of this paper is to give a method for a computation of a 
non-perturbative running in the operators mixing case. 
Furthermore the operators that we consider here already 
have an important phenomenological relevance since 
they allow for the computation of standard model and 
beyond the standard model neutral kaon mixing matrix 
elements and of iiT — tttt amplitude in the A/ = 3/2 
channel. 



Numerical application 

The RBC-UKQCD collaboration has recently performed 
a computation of if — > tttt decay amplitudes (in both 
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FIG. 2: Renormalization factor Zi\/Z\ — Z^^, obtained at 
two different values of the lattice spacing. Except for a few 
points, the error bars are smaller than the symbols. 



isopsin channels) with a two-pion final state [25]. The 
physical matrix elements are computed from the Eu- 
clidean ones by using the Lellouch-Liischer formula 1261 . 
The results are very promising but, since it was the first 
computation of its kind, an unphysical pion mass of 
rriT^ ^ 422 MeV was used. In order to simulate the phys- 
ical kinematics the collaboration is currently repeating 
the computation on a much larger volume - of spatial 
extent Lq ~ 4.6 fm - with a nearly physical pion mass of 
•niT, ~ 140 MeV. Some promising preliminary results of 
the if — >■ TTTT matrix elements in the A/ = 3/2 channel 
have been reported in |27|. The gauge action used is a 
modification of the Iwasaki gauge action following the 
lines of [28, ,2£, , that we call the "dislocation suppressing 
determinant ratio" . On the same lattice, a computation 
of the matrix element (if lOyy^tj^^^jif is also on the 
way. With a lattice spacing ao ^ 0.14fm one might 
worry about the size of the discretisation effects for a 
momentum of two or three GeV (indeed one might doubt 
the existence of the Rome-Southampton window on this 
lattice). The way out is to follow the strategy explained 
in the previous section with n = 1 (two different physical 
volumes). Firstly, in the volume Lq we compute the 
bare matrix elements and the renormalization factors 
at a low energy /io, where we use /io ~ 1.5 GeV. Then 
the continuum limit of the step scaling matrix cr(/ii, /io), 
with /il €E [/io, 3 GeV] is obtained from another set 
of simulations. There we use Li 2.7 fm and two 
finer lattice spacing ai ~ 0.086 fm and 0.114 fm (more 
details about these simulations can be found in [30 1 131 ] ). 
Finally we plan to combine the two results in order 
to compute the renormalized matrix element at 2 or 
3GeV and then apply the perturbative matching to 
MS. As mentioned above the two different volumes 
use two different gauge actions but the results can be 
combined together because we extrapolate the step 
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scaling matrix to the continuum. Although in principle 
we would like to have a third lattice spacing in order to 
have a better control on the continuum extrapolations 
in the volume Li, in practice the discretisation effects 
on the step scaling matrix elements appear to be small. 
Obviously this method can be applied to all sort of 
different quantities that one would like to extract from 
our large lattice Lq. In the near future, when the next 
generation of supercomputers will be available, we plan 
to add a finer lattice on the large volume Lq. 

In the rest of this section we present our results for the 
step scaling matrices a{fio,fi) with fi varying in a range 
[/X0i~ 3GeV]. The computation of the renormalization 
factors on the fine lattice is done using the same setup 
as in [131 [32], to which we refer the reader for a more 
detailed explanation. The computation of the Z factors 
has been already presented in [32j . One of the main fea- 
tures of our computation is the use of the Domain- Wall 
fermion 33 35J which exhibits an almost exact chiral- 
flavour symmetry. As a consequence the renormalization 
pattern is the same as in the continuum (up to some 
numerically irrelevant lattice artefacts). Some other in- 
teresting aspects of our calculation are the use of the 
volume sources [36] giving us a very good statistical pre- 
cision, non-exceptional kinematics [TU] [H] to suppress 
the Goldstone pole contributions: here we use a scheme 
which is called (7^, 7^)-scheme in [T3]. In figure [2] we 
show the Z factor of the (27, 1) operator, normalized by 
Zp^ and extrapolated to the chiral limit, obtained on the 
two different lattices at the simulated momenta (covering 
a range from ~ 1.1 GeV to ^ 3.5 GeV). The correspond- 
ing step scaling function is computed according to the 
definition eq Q. In figure [3] (left) we plot the result at fi- 
nite lattice spacing, together with the continuum extrap- 
olation. In order to determine the vertex function at any 
given momentum for each ensemble we fit our data as a 
function of and interpolate. As explained earlier in the 
text, because we are using twisted boundary conditions 
that keep the momentum orientation fixed with respect 
to the lattice axes this is a smooth interpolation. Hav- 
ing obtained in this way the vertex function on a number 
of lattice spacings at a fixed physical p"^ we extrapolate 
this to the continuum limit using a constant plus 0{a^) 
ansatz for each value of the momentum. This ansatz is 
justified in our approach because as we pick a fixed mo- 
mentum orientation with respect to the lattice axes, 0(4) 
breaking lattice artefacts are well parametrised [1^. As 
one can see on the plot, the lattice spacing dependence 
is very well under control. In ,32] we have shown that 
the renormalization pattern is the same as in the contin- 
uum, the chirally forbidden renormalization factors being 
zero within error. The chiral extrapolation is done using 
three different quark masses and we found a very mild 
quark mass dependence for all our quantities. In our 
setup the sea quark mass of the strange is fixed (to its 
physical value), whereas the valence sea quarks masses 
are equal to the sea light quark masses and then extrap- 



olated to zero. As a consequence our results are affected 
by a small systematic error, which was evaluated in |13j 
for the (27, 1) operator. We have checked that taking 
the chiral extrapolation of the ratio of Z gives the same 
result as taking the ratio of the chiral extrapolation of 
Z. In figure |3] (right) we compare our non-perturbative 
result with the next-to-leading order (NLO) running [13] . 
We note that the running is quite small (~ 4% between 
fi ^ 1.5 GeV and ~ 3.5 GeV), while NLO perturbation 
theory predicts ~ 7%. With our very small error bars 
such a the difference is clearly visible. In figure |4j we 
plot our results for the two (continuum) two-by-two step 
scaling matrices cr(/i,/io) computed in the renormaliza- 
tion basis. By definition, at the matching point /x = /xq 
the matrices are equal to the identity. In figure [5] we 
compare our results for the electroweak penguins in the 
As — 1 basis to the next-to-leading order (NLO) run- 
ning recently computed in |37j . In Figure [6] and [7] we 
set fxo — 3GeV and plot the non-perturbative running 
divided by the NLO prediction. In general we find that 
when NLO perturbation theory is available the results 
agree qualitatively with the non-perturbative ones. It is 
interesting to note that for both the (8, 8) and the (6, 6), 
we found that one off-diagonal matrix element has a very 
small non perturbative running, and that the correspond- 
ing one loop anomalous dimension is either zero (for (732) 
or a very small number (for (745). 

Conclusion 

We have presented a general method to compute the 
non-perturbative continuum running in the operator 
mixing case. In particular this method allows the use 
of an RI-MOM type scheme on a rather coarse lattice. 
We have computed this running between ^ 1.5 GeV and 
~ 3.5 GeV in the case of four-quark operators which 
occur in neutral kaon mixing (including the BSM ones) 
and AI = 3/2 K ^ nn decays. Although the strategy 
we have presented is very general we have shown why it 
is important for the light quark physics program of the 
RBC-UKQCD collaboration and in particular for a full 
computation of K — > tttt amplitudes. 
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FIG. 3: Left: step scaling function (non-perturbative running) of the (27, 1) operator at finite lattice spacings and in the 
continuum in the lattice SMOM (7^, 7fi)-scheme . The energy scale ^0 is fixed to fio ~ 1.5 GeV, and ^ varies in the range 
[1. GeV, ~ 3.5 GeV]. The errors bars are smaller than the symbols. Right: the same quantity is compared to perturbation theory 
(dashed black curve : one loop, solid black curve: two loops, solid blue curve: continuum extrapolation of non-perturbative 
running with its error) in the range [^0, ~ 3.5 GeV]. 
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FIG. 4: Step scaling matrix cr(/i, /iq) given in the As — 2 renormalization basis. On the left panel we show the matrix elements 
of the (8,8) operators and on the right the (6,6) operators. Like in figure [s] the energy scale fio is fixed to fio ~ 1.5 GeV, and 
fi varies in the range [fio, ~ 3.5 GeV]. 



APPENDIX 



1. Denoting by the superscript "+" the parity even 
component, the expUcit relation between the two As = 2 
basis is 



(6,6) 




- Q 



As=2 
3 : 

2Q2 " 



(27,1) [Of^=2]- 
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2. To define the AJ = 3/2 part of the As = 1 operators, we follow the conventions of [ID] 

Qi = (5^7^(1 - 75)^0) [(u/STa'CI ~ 75)^/3) - ('?^)7p(l -75)^/3))] + (sa7A<(l - 75)mq)('«/37m(1 - 75)^/3)) (A.l) 

Q'l = (Sa7A.(l - 75)^0) [(",97m(1 +75)'"/3) - (s^)7m(1 +75)S;9))] + (Sa7^(l -75)Wa)(u/}7^(l +75)d^)) (A.2) 

Qs = (sa7M(l - 75)^/3) [(w^7m(1 + 75)^0) - (s^7^(l + 75)sq))] + (sq7^(1 - 75)u^)(u^7^(l + 75)da)) (A. 3) 
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FIG. 5: Step scaling matrix a{iJ.,iJ,o) of the (8,8) electroweak penguins in the As = 1 renormalization basis. For each matrix 
element we compare to perturbation theory (dashed black curve : one loop, solid black curve: two loops, solid coloured curve: 
continuum extrapolation of non-perturbative running with its error). 



8 



1.01 



a. 



0.99 



3 

^0.98 
3- 0.97 
b 0.96 
0.95 




5 6 7 



10 



FIG. 6: Step scaling function (Tii(/^,/io) of the (27, 1) operator divided by the corresponding perturbative running Und-i, no) 
(dashed line: one loop, solid coloured curve: two loops). For ease of comparison with perturbation theory and in contradiction 
to the previous plots, is fixed at the conventional scale of 3 GeV and fj. varies in the range [1.5 GeV, /.to]. 
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FIG. 7: Same as figure [s] for the Step scaling matrix a{fj., fj,o) of the (8,8) electroweak penguins. 
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